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ABSTRACT 

We present temperature and metallicity maps of the Perseus cluster core obtained 
with the Chandra X-ray Observatory. We find an overall temperature rise from ~ 3.0 
keV in the core to ~ 5.5 keV at 120 kpc and a metallicity profile that rises slowly from 
~ 0.5 solar to ~ 0.6 solar inside 60 kpc, but drops to ~ 0.4 solar at 120 kpc. Spatially 
resolved spectroscopy in small cells shows that the temperature distribution in the 
Perseus cluster is not symmetrical. There is a wealth of structure in the temperature 
map on scales of ~ 10 arcsec (5.2 kpc) showing swirliness and a temperature rise 
that coincides with a sudden surface brightness drop in the X-ray image. We obtain 
a metallicity map of the Perseus cluster core and find that the spectra extracted 
from the two central X-ray holes as well as the western X-ray hole are best-fit by gas 
with higher temperature and higher metallicity than is found in the surroundings of 
the holes. A spectral deprojection analysis suggests, however, that this is due to a 
projection effect; for the northern X-ray hole we find tight limits on the presence of 
an isothermal component in the X-ray hole, ruling out volume-filling X-ray gas with 
temperatures below 11 keV at 3er. 

Key words: galaxies: clusters: general - galaxies: clusters: individual: Abell 426 - 
cooling flows - intergalactic medium - X-rays: galaxies 



1 INTRODUCTION 

The Perseus cluster, Abell 426, at a redshift of z = 0.0183 or 
distance about 100 Mpc is the the nearest high lu minosity 
cluster with a high central su r face brightne ss (e.g. 



not distinguishable as sharp features on a 3-7 keV image 
and are there fore not shock features, contrary to the early 
prediction of Heinz, Reynolds & Begelman (1998), and are 



et al. 1981; Allen et al. 2001 ; Fabian 1994). This makes it 
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the brightest cluster in the X-ray sky. 

The first Chandra subarcsecond-resolution X-ray im- 
ages of the cluster core aroun d the central domin ant galaxy 
NGC 1275 were presented by |Fabian et al. (2000|, F00). Us- 



ing X-ray colours F00 showed that the temperature of the 
gas decreases from about 6.5 keV to 3 keV inward to the 
central galaxy NGC f 275, which is surroun ded by a spe ctac- 
ular low-ionization, emission line nebula ( Lynds f 97C , see 
also Conselice, Gallagher, & Wyse 2001). 



not expanding supersonically. The simplest interpretation of 
the low surface brightness is that they are devoid of X-ray 
gas and have pressure support from cosmic rays and mag- 
netic fields. However they may contain some hotter gas at 
the virial temperature, or even above, of the cluster with the 
radio plasma having a low filling factor. In this scenario the 
rims consist of cooler gas which has been swept aside. The 
two outer holes, one of which was previously known, were 
identified w ith recently-found sp urs of low-frequency radio 
emission by Blundell et al. (200C ). 



The nucleus powers the radio source 3C84 (Pedlar et 
al f990) which has structures on various scales. The 0.5 
arcmin-sized radio lobes coincide with holes in the soft X- 
ray emission ( Bohringer et al. 1993) ; McNamara et al. 1996 
F00). Chandra has revealed that the holes are not due to ab- 
sorption and have X-ray bright rims, which are cooler than 
the surrounding gas (F00). It was found that the rims are 
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The core of the Perseus cluster offers the rare oppor- 
tunity to study the radio source/intracluster gas interaction 
on arcsecond scales. In the present study we have carried out 
a spectral analysis of the Chandra images which provides us 
detailed information about the cluster physics. In Sect. |we 
describe the observations including the full 29.0 ks Chan- 
dra image of the Perseus cluster core. In Sect. [] we describe 
our method of spectral analysis and show temperature and 
metallicity profiles, as well as maps. In Sect. ^|we carry out 
a spectral deprojection analysis of the cluster and search 
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for the presence of X-ray gas in the northern X-ray hole. In 
Sect. ^ we conclude with a summary and a discussion. 

We use Ho = 50kms _1 Mpc -1 and go = § throughout. 
Unless otherwise stated, quoted error bars are 1 a (68.3% 
confidence). 



2 OBSERVATIONS 

Chandra observed the Perseus cluster on 1999 September 20, 

1999 November 28 and 2000 January 31 for 5.3 ks, 9ks and 
24.7 ks, respectively. The early observation was carried out 
with the ACIS-I detector. The energy resolution of this ob- 
servation has been seriously affected by the radiation dam- 
age suffered by the front-illuminated ACIS chips, so that 
we only use it to produce images. The two later observa- 
tions were carried out with the ACIS-S detector. The target 
was centred very close to the middle of the back-illuminated 
ACIS chip S3 (ACIS-S3). The detector temperature during 
these observations was -110 C. Inspection of the light curves 
of the observations showed that the later, longer, observa- 
tion had been affected by flares, so that we had to clean 
the light curve. This reduced the exposure time usable for 
spectroscopic analysis to 14.9 ks. 

In Fig. [l] the combined 29.0 ks image of all three Chan- 
dra exposures of the Perseus cluster core is shown. The image 
was constructed by concatenating the three event files using 
the COMBlNE_OBSlD script available on the Chandra X-ray 
Observatory Centre (CXC) web pages (www.cxc.edu). 

The image was binned into (1.97 arcsec) 2 pixels. The 
well-known features of the inner X-ray holes and the outer 
X-ray hole (north-west) can be seen. In this stacked image a 
streak of emi ssion can be seen r unning through the southern 
hole (see also Fabian ct al. 2002, hereafter F02). This feature 
was not obvious in the image studied by F00 due to the dead 
column feature that runs straight across it in the January 

2000 image. The southern outer X-ray hole first discussed 
by F00 and another similar cavity immediately to the west 
of it are also visible. 

Between the western outer hole and the southern outer 
hole a curved, elongated feature is apparent with a smaller 
count rate than in the surrounding area. It appears to be a 
structure that connects the two outer holes. It is reminiscent 
of the radio structure apparently connecting the two outer 
radio lobes of the central radio source 3C 338 in Abell 2199 



(Giovannini et al. 1996 



The appearance of swirliness (F00, see also Churazov et 



al. 2000) in the outer regions of the image is created by a 



distinct surface brightness drop that broadly coincides with 
a rise in the X-ray temperature (with increasing radius), as 
we will show in the next section. 



3 SPECTRAL ANALYSIS 
3.1 Calibration 

We extract spectra from regions of the Chandra data set us- 
ing the CI AO software package distributed by the CXC. We 
only work with photons detected in the ACIS-S3 chip, which 
is also known as ACIS chip 7. We also restrict ourselves to 



the 0.5-7.0 keV band in order to be certain of a reliable en- 
ergy calibration. Response matrix files (RMF) and ancil- 
lary response files (ARF) were produced using the calcrmf 
and calcarf programs by Alexey Vikhlinin available from the 
CXC web pages. 

Using the task GRPPHA from the FTOOLS software 
package provided by the NASA High Energy Astrophysics 
Science Archive Research Center, the extracted spectra are 
grouped so that we have at least 20 counts in a bin corre- 
sponding to a certain range of energy channels. This allows 
us to use Poisson error bars for spectral fitting. 

We have have generated background spectra using 
Markevitch's method and program that is also available from 
the CXC web site. We extract the background data from the 
same regions on the chip as the spectra extracted from the 
data sets. In the case of the Perseus cluster core the back- 
ground in the ACIS-S3 chip in the 0.5-7.0 keV band (< 1 
ctss -1 ) is much smaller than the signal from the cluster 
(83.4 ctss" 1 ). 

When fitting extracted spectra, we fit the November 
and the January data set in parallel with the same spectral 
model, allowing only for a different normalization of the two 
spectra. The different normalizations are necessary because 
the exposures are affected by the ACIS-S3 missing column 
feature in two different orientations. The feature is due to 5 
contiguous partially dead columns in the centre of the chip 
that cause missing counts in an elongated region several ar- 
cmin long and approximately 17 arcsec wide (due to the 
spacecraft dither) . This procedure affected only regions cov- 
ering the dead columns, in other regions the normalizations 
were consistent with each other. 



3.2 Selection of regions 

In order to map the spectral properties of the cluster we 
divide the cluster into cells with a constant number of cluster 
emission count s. S ince the background is not important for 
Perseus (Sect. 3.1) this is very similar to having an equal 
number of counts in each cell. 

We determine the cell shapes by iteratively dividing an- 
nular segments with inner and outer radii n, and 12 and 
bounding azimuthal angles ai,«2 into regions with equal 
number of counts (corrected for the background contribution 
as determined from Markevitch's background files). Begin- 
ning with a full circle with radius r (ri = 0, r-z = r, ol\ — 0, 
«2 = 2 7r) the iterations are either exclusively performed on 
the radius, yielding a sequence of annuli, or interchangingly 
on radius and azimuth, yielding a more complex web of cells. 

We carry out the analysis on three scales: (1) Simple an- 
nuli with ~ 35000 counts each, (2) cells with ~ 18000 counts 
and (3) cells with ~ 500 counts each. The cell divisions were 
determined using the 14.9 ks exposure taken on 1999 Novem- 
ber 28. The mentioned count values in this and the following 
sections refer to this exposure. 



3.3 Annular spectra 

We begin the spectral analysis by studying azimuthally av- 
eraged spectra from annuli with a large (~ 35000) number 
of counts. This will give an overview of the radial variation 
of the ambient X-ray gas temperature in Perseus. 
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Figure 1. Combined 29.0 ks Chandra image of the Perseus cluster core in the 0.5 to 7 keV band. The image size is 10.1 arcmin X 7.9 
arcmin (312.5 kpc X 244.4 kpc). The linear boundaries visible at the edges of the frame are due to the three different orientations of the 
detector in the three added exposures. 



We use XSPEC ( |Arnaud 1996] ) to fit the X-ray spectra. 
Errors are determined with the XSPEC tasks ERROR and 
STEPPAR. We model the X-ray emission from the intraclus- 
ter gas using the APEC plasma emission code (version 1.10). 
It is available on the web at hea-www.harvard.edu/APEC 
and is also included in the latest XSPEC versions. We use 
this n ewer code instead of the traditional MEK AL plasma 
code (Kaastra & Mewe 1993; Liedahl et al. 1995) because it 



provides a significant improvement of the spectral fits. For 
example, when fitting a spectrum extracted from an annu- 
lus with inner radius 23 arcsec and outer radius 32 arcsec 
around the nucleus (542 degrees of freedom) we obtain a 
X 2 = 625.7 for the APEC model instead of X 2 = 662.0 
for MEKAL. The change in best-fit parameters, however, is 
small; the best-fit absorption and temperature are identi- 
cal for both plasma codes. The metallicity changes slighly, 
in this particular case from 0.42±0.03 solar for MEKAL to 
0.46T0.03 solar for APEC. A large part of the improvement 
is due to a better representation of the blue wing of the Fe 
L complex. 



In order to model the emission from a single tempera- 
ture we fit the spectra with the following model: 



MODELi = PHABS (jV H ) x APEC (T, Z, K) 



(1) 



PHABS is the photoelectric absorp tion model by 



Balucinska-Church & McCammon (1992). The free param- 
eters equivalent hydrogen column density Ah, temperature 
T, ave rage metallicity Z (relativ e to the solar values accord- 
ing to Anders & Grevesse 1989 ) and spectral normalization 
K are given in brackets. In Figs. ^, ^ and ^ the results of the 
spectral fits for the temperature, the metallicity and for the 
required photoelectric absorption are shown. The fits have 
been done both with (dotted error bars) and without (solid 
error bars) fixing the photoelectric absorption. Except for 
the innermost annulus, the absorption is consistent with 
being constant (Fig. ^). For the fits with fixed absorption 
we fixed it at the median of all fitted equivalent hydrogen 
column densities 7Vh = 0.143 x 10 22 cm~ 2 because the 



Galactic value (Dickey & Lockman 199C) value from Hi 
studies (0.149 x 10 21 cm" 2 ) appears to be slightly too high 
(see Fig. §). In Tab. [j] we list the details of the fits with 
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free absorption. It can be seen from Figs. ^| and ^| that the 
best-fit values are very similar for models with or without 
fixed absorption. 

The temperature profile (Fig. ^ reveals a rapid temper- 
ature decline from more than 5 keV to about 3 keV in the 
rims around the inner X-ray holes. Due to the large number 
of counts, the error bars are very small. There is a sudden 
jump in the temperature at a radius of 80kpc, which coin- 
cides with the sudden brightness drop seen to the north of 
the cluster centre (Fig. Q). 

Note also that the temperature rises again in the core. 
In this innermost ring, the nucleus and the inner 2.5 kpc 
(4.8 arcsec) around it have been cut out. It can be seen from 
Tab. [l] that the spectral fit in this ring is not as good as 
the fits in most of the other rings. This region encompasses 
most of the northern and about half of the southern X-ray 
hole. It shows signs of excess a bsorption (Fig. ^ ) due to the 
infallin g high-velocity system ( Unger et al. 199o| ; |de Young 
et al. 1973; Briggs et al. 1982) that was found in absorption 
by Chandra (F00). A two-temperature model improves the 
fit and is favoured by an F-test. The two-temperature fit 
yields a 4.6 keV gas phase and a 1.9 keV gas phase, where 
the normalization of the high-temperature phase is twice 
of the low-temperature phase. It only improves the reduced 
Xv 2 to 1.34, however, so that this can also not be viewed 
as an entirely acceptable fit. In general we find that two- 
temperature fits do not provide a better description of the 
annular spectra. Overall there is no convincing evidence for 
multiphase gas in the Chandra spectra of the Perseus cluster. 

In the metallicity profile (Fig. ^) a gradient from Z ~ 
0.6 solar at a radius of 60 kpc to Z ~ 0.4 solar at a radius 
of 120 kpc is seen. There may be a peak of the metallicity 
profile at ~ 50 — 60 kpc, but the metallicities drop only 
slightly inwards from 60 kpc towards a central metallicity of 
Z ~ 0.5 solar. Due to the large number of counts in the 
annuli, the relative errors are normally less than 10 per cent 
so that the trend between 60 and 120 kpc is detected with 
high significance. 



3.4 Spectra in large cells 

Following on to the azimuthally averaged spectra we next 
present the results from cells with ~ 18000 counts. In the 
left panel of Fig. |E| the definition of the cells and the inten- 
sity map of the cluster with this resolution are shown. The 
contours of the Chandra image are overlaid on top. The mid- 
dle and right panels are maps of the temperature and the 
average metallicity determined using the single temperature 
model in eq. (|l]) with free Galactic absorption. The statisti- 
cal uncertainty of the temperatures is less than 0.1 keV, the 
uncertainty of the metallicities is less than 10 per cent. 

It was shown by F00 using X-ray colours that the tem- 
perature distribution in the region is not circularly symmet- 
ric. The middle panel in Fig. ^illustrates the overall temper- 
ature distribution obtained with the APEC plasma model 
fits in a region of ~ 125 kpc radius around the nucleus. The 
large asymmetric region in the centre at a temperature of 
3 to 4keV is surrounded by hotter gas with temperatures 
up to 6keV. The cool region is extended towards the north 
and the west, whereas in the south-east hotter gas is found 
at smaller radii. A comparison with the intensity contours 
in the left panel shows that the outer intensity contour has 
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Figure 2. Temperatures from single temperature fits to annular 
spectra between 2.5 and 121.8 kpc (see also Figs. ^ and The 
temperatures with solid error bars have been determined with- 
out fixing the galactic absorption (see Fig. ^). For the dotted 
error bars the absorption was fixed to the median of the fitted 
absorption values, Nh=0.143x10 22 cm -2 . Note the jump in the 
temperature profile at a radius of 80 kpc. 
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Figure 3. Metallicity profile from single temperature fits to annu- 
lar spectra between 2.5 and 121.8 kpc (see also Figs. |and|). The 
metallicities with solid error bars have been determined without 
fixing the galactic absorption (Fig. ^). For the dotted error bars 
the absorption was fixed to the median of the fitted absorption 
values, N H =0.143xl0 22 cm -2 . 
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Table 1. Single temperature fits to azimuthally averaged spectra. The listed parameters are inner and outer radius (ri, r 2 ), temperature 
(column 3), metallicity (column 4) and equivalent absorbing hydrogen column density (column 5). Describing the goodness of fit we also 
list the x 2 value, the number of degrees of freedom (DOF) and the reduced \v 2 (columns 6 to 8). 
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a similar shape to the isotherms at that radius; the inward 
(towards smaller radii) bend of the outer intensity contour 
to the south east of the cluster centre coincides with a high- 
temperature region. The spectra extracted from the regions 
corresponding to the inner X-ray holes and the outer hole 
(see the lighter cell the left panel of Fig. |^) in the west are 
best-fit by gas at a higher temperature than is found in their 
immediate surroundings. The same may be true for the less 
obvious southern X-ray hole discovered by F00 (see Fig. ^), 
but here the distinction from the hotter gas that intrudes 
from the south-east is not clear. 

In the metallicity map (Fig. ^, right panel) small varia- 
tions of the metallicities are apparent between cells at sim- 
ilar radii. A high-metallicity ring (~ 0.6 solar) around the 
cluster centre is prominently seen at a radius of ~ 60 kpc 
where the peak in the azimuthally averaged metallicity dis- 
tribution in Fig. |^ was found. In addition to the higher tem- 
perature, the spectra extracted from the inner X-ray holes 
and the western hole are also best-fit by gas with a higher 
metallicity than is found in their immediate surroundings. 

We note that a comparison with the galaxy distribu- 
tion in an optical Digitial Sky Survey image of this region 
does not yield any obvious correlation. Due to the relatively 
coarse scale of the metallicity map, and due to projection 



effects, however, the question of correlations between the 
metal distribution and the galaxy distribution cannot be 
answered by this data set. Nevertheless, it is exciting that 
the Chandra X-ray observations begin to approach answer- 
ing such questions. 

3.5 Spectra in small cells 

In this section we finally use two grids of cells: firstly we 
determine the amount of absorption and the metallicities in 
cells with ~ 9000 counts each using the model in eq. We 
then subdivide these cells further until we obtain a grid with 
~ 500 counts in each cell. In order to obtain cells that are 
not too elongated, we limit the ratio 

(a 2 - ai) x (n + r 2 )/2 



r 2 - n 



(2) 



of cell width and depth (see Sect. 3.2) to be greater than 
1/5 and smaller than 5. The cells have a typical size (circle 
section at half radius) of 10 arcsec (5.2 kpc). In these fine 
cells we determine the temperature with a fixed metallicity 
at the value determined in the larger cells. This procedure 
yields the finely meshed temperature map shown in Fig. [j| 
The relative errors of the temperatures in this plot are of 
the order of 10 — 15 per cent. 
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Figure 4. Column densities from single temperature fits to annu- 
lar spectra between 2.5 and 12L8 kpcwithout fixing the absorbing 
column density (see also Figs. Hand H). The dashed line indicates 
the median of N H =0.143xl0 22 cm-^. 



In order to obtain stable solutions for the absorption, we 
fix it to the median value Nh = 0.143 x 10 22 cm" 2 (Fig. f| if 
it is smaller than A^h = 0.16 x 10 22 cm -2 , or otherwise fix it 
at the value determined in the first step; in total, there were 
only two cells where the local absorbing column densities 
of 7Vh = 0.20 x 10 22 cm" 2 and N H = 0.21 x 10 22 cm" 2 , 
respectively, had to be used. These two cells correspond to 
the region where the infalling high- velocity system is seen in 
absorption in the X-ray colour image by F00. 

When looking at the goodness of the individual fits 
we noticed in a small number of cases that the assumed 
(fixed) metallicity appeared to underestimate the metallic- 
ity in these cells. In the cell with n = 43.39 kpc, V2 = 
46.43 kpc, Qi = 346° and oti = 351° (angles given in de- 
grees, counterclockwise from north), for example, we found 
a reduced x^ 2 value of 1.3 when the metallicity was fixed 



Figure 6. Temperature map of the perseus cluster core with 
500 counts/cell. T he c ounts were determined using the 14.9 ks 
exposure (see Sect. 3.2). The coordinates are J2000. 



at Z — 0.434 solar. The reduced Xv drops to 1.1 for the 
best fit Z — 2.87j^'®4 solar. The constraint is entirely due 
to the iron L complex. This also leads to an underestimate 
of the temperature: 3.8±o gg keV compared 4.57± ;f| keV 
if the metallicity is left as a free parameter. Although the 
evidence is not statistically compelling, this may be due to 
small-scale inhomogeneities of the metallicity. However, cur- 
rently the count rates on small scales are too low and ac- 
cordingly the uncertainties too large to reliably determine 
metallicities on these scales. 

In the temperature map the overall temperature distri- 
bution from Fig. ^ (middle panel) is again visible. In addi- 
tion, the map includes fine detail showing prominently that 
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the temperature is higher in the regions of the two inner 
and the north-western X-ray holes than in their surround- 
ings. Extending from the bright rim to the south-east of the 
central X-ray holes the low-temperature (2 to 3 keV) gas ex- 
tends westwards and spirals around counter-clockwise. This 
spiral shape of the colder gas corresponds to the swirly ap- 
pearance of the 0.5 to 7 keV emission shown in Fig. jl]. The 
shape of lines of constant temperature at radii ~ 80 kpc (2.6 
arcmin) discussed in Sect. ^4 are also seen here to be simi- 
lar to the shape of the surface brightness drop in the X-ray 
emission. The general spiral appearance of the cooler X-ray 
gas in the core was also noted by F00 in their X-ray colour 
image; the shape of the temperature distribution and the 
overall magnitude of the temperatures in the cluster core as 
determined with X-ray colours by F00 are confirmed by our 
individual spectral fits. 



4 SPECTRAL DEPROJECTION 

4.1 Azimuthally averaged profiles 

In the previous sections we have treated the cluster as a two- 
dimensional object in order to learn about the overall spec- 
tral properties. Since complicated projection effects must be 
present due to the different locations of, for example, bub- 
bles, brightness drops or metals in the cluster, we attempt 
in this section to unravel some of the projection effects that 
are present. However, this can only be done by assuming a 
specific symmetry. 

We deproject the observed spectra, in particular the 
temperature profile (Fig. ||) and the metallicity profile 
(Fig. ^|), by assuming spherical symmetry and using the 
PROJCT model available in XSPEC combined with our stan- 
dard APEC model (eq. |l]) . In order to reduce the number of 
parameters, only eight annuli with ~ 144, 000 counts each 
were used (as measured in the January 2000 image). 

The model assumes eight spherical, concentric shells 
filled with gas with independent temperatures and metal- 
licities. The projct model adds the APEC spectral contri- 
butions from the X-ray emission in the shells in order to 
model each of the eight observed annular spectra. All eight 
spectra are fitted at the same time^]. Consistent with Fig. 
we find no excess absorption except for the innermost region. 
With the current model and due to the complicated shape 
of the absorbing area (F00), however, it is not possible to 
determine in which shell the absorbing material is located. 

The final deprojected temperature and metallicity pro- 
files are shown in Fig. [5]. We plot the deprojected values with 
solid error bars, while the temperatures and metallicities de- 
termined from the projected spectra are plotted with dotted 
error bars. Starting from the outside, the temperature pro- 
file drops very quickly in the first bin and continues to drop 
monotonically, but with a shallower slope afterwards. Also 
starting from the outside, the metallicity gradient appears 
to be steeper down to ~ 70 kpc than it was for the pro- 
jected spectra, but is very similar, albeit with larger error 



The PROJCT model only deprojects a single data set at a time. 
We thus used only the longest 14.9 ks data set from January 2000 
for the deprojection analysis reported in this section. 
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Figure 8. Electron pressure profile as determined using a spec- 
tral deprojection of the Perseus cluster core. The sudden rise in 
pressure in the outermost bin is not real, but an edge effect due 
to our data analysis (see text). 



bars, within 70 kpc. Inside a radius of ~ 70 kpc, we find the 
metallicity approximately constant at Z ~ 0.5 solar. 

The X-ray data also allow us to determine the emis- 
sion measure EM = x V (via the normalization constant 
K in eq. [j]), where n c is the electron density and V the 
volume of the gas. Assuming that isothermal X-ray gas at 
the deprojected temperature T fills the spherical shells, we 
can determine the pressure p of the electrons, p — n e kT, 
as shown in Fig. ^. Note that the sudden rise in pressure 
in the last bin is an edge effect due to our assumption that 
there is no cluster emission outside this radius. Similarly, the 
steep deprojected temperature drop in the outermost bin of 
Fig. ^ is an artifact since the cluster extends significantly 
beyond the outermost radius accessible to the ACIS-S3 de- 
tector used. This is also apparent from the fact that the 
deprojected temperature value is best-fit higher than the 
projected value. The temperature, metallicity and emission 
measure values further in, however, will not be affected by 
this due to the steep brightness increase (see Fig. ^, left 
panel) towards the nucleus. 



4.2 The X-ray emission from the X-ray holes 

In the projected temperature maps we found that the spec- 
tra extracted from the two central X-ray holes as well as the 
western X-ray hole are best-fit by hotter and more metal 
rich gas than their immediate surroundings. It is possible 
that this is due to a projection effect, whereby we would be 
seeing the hot gas above and below an empty cavity. 

In order to investigate the properties of the X-ray gas 
that fills the region of the holes, we have carried a local 
deprojection analysis. We have spectrally deproj ecte d the 
emission emitted by 9 annular segments (Sect, 
azimuthal boundaries ai = 350 degrees and 02 



3J) with 
28 de- 
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Figure 7. Deprojected temperature (left panel) and metallicity profiles (right panel) of the Perseus cluster core (solid lines). Also plotted 
with dotted lines are the corresponding projected quantities that are shown in Figs. ^ and |i] with higher resolution. 



grees and 14000 counts per segment (as measured in the 
January 2000 image) in the 0.5 to 7 keV band. The in- 
nermost segment occupies the region of the northern inner 
X-ray hole with radii between 5 and 27 arcsec of the nu- 
cleus. The analysis was analogous to the full-circle analysis 
in the last section, except that the absorption was forced 
to be the same for all segments and the metallicity in the 
inner hole was fixed at Z = 0.5. The best-fit absorption was 
N H = 1.40 x 10 21 cm -2 , the reduced \ 2 °i the best fit was 
xl = 1-14. 

Interestingly, we find that almost the entire projected 
emission in the innermost segment can be explained by the 
emission from the shells above. In Fig. ^| we show the al- 
lowed 68% (1 a), 95.4%(2 a) and 99.73%(3 a) confidence 
contours in the temperature/electron density plane of the 
parameters for an APEC component in the innermost seg- 
ment. It can be seen that only a very low density plasma 
is allowed by the analysis. Since the emission measure con- 
strained by the observed X-ray data is proportional to the 
volume of the gas, the allowed electron densities depend on 
the filling factor / of the X-ray gas in the hole. We thus 
show the lines of pressure equilibrium with the second vol- 
ume segment for three different filling factors of the X-ray 
gas (the electron pressure in the deprojected second segment 
is p — n c T — 0.17keVcm -3 , slightly higher than the shell- 
averaged pressures given in Fig. ^|). These data rule out at 
3cr the presence of X-ray gas at a temperature lower than 
11 keV filling the entire X-ray hole. At the same confidence 
level we can rule out intracluster gas at the virial tempera- 
ture of the cluster (~6.5 keV), filling one third of the X-ray 
hole. 




6xl0" 3 8xl0- 3 10" 2 1.2xl0" 2 1.4xl0" 2 1.6x10-= 
■v/f x electron density (cm -3 ) 



Figure 9. 68% (1 cr), 95.4%(2 a) and 99.73%(3 a) confidence con- 
tours (thick lines) on the presence of an isothermal APEC com- 
ponent in the northern inner X-ray hole. The hyperbolas (thin 
lines) denote the required parameter pairs for pressure equilib- 
rium with the X-ray gas at larger radii, assuming filling factors 
f = jfi, f = 1 and / = 1 of the X-ray gas. 



5 DISCUSSION 

We have presented in this paper the results from a detailed 
temperature and metallicity mapping of the Perseus cluster 
core using Chandra observations. 

We work with the complete useful ACIS-S data set com- 
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prising a total of 23.9 ks. Beginning by averaging over large 
annuli around the nucleus, we find that the temperature av- 



aclustcr gas has recently become the subject of much work 



eraged n such annuli rises smoothly from ~ 3.0 keV to ~ 5.5 



hurazov et al. 2001 



& Balogh 2001 



keV at 120 kpc. A small kink is found, however, in this pro- 
file at ~ 80 kpc, which corresponds to a surface brightness 
drop seen in the X-ray image. Using solar abundances ac- 
cording to Anders fc Grevesse (1989| ), the metallicity rises 
from Z ~ 0.4 solar with decreasing radius to a maximum 
Z ~ 0.6 solar at ~ 60 kpc. There may be a peak of the 
metallicity profile at ~ 50 — 60 kp c, which is reminiscent o f 
such peaks in other clusters (e.g., Sanders fc Fabian 2001 ). 
However, inside 60 kpc the metallicity varies little, with a 
central metallicity of Z ~ 0.5 solar. We have also carried 
out a spectral deprojection analysis of the radially averaged 
profile, which yields a very similar picture. 

Spatially resolved spectroscopy in small cells shows that 
the temperature distribution in the Perseus cluster is not 
symmetrical. In fact, the distribution of cold (2 — 3 keV) 
gas in the centre appears to spiral outward and corresponds 
to the swirly appearance of the X-ray emission (see Figs, jj] 
and ^ . This may suggest the presence of angular momentum 
of the intracluster gas (F00). T here are, however, other pos- 
sibilities, such as the model by |Churazov et al. (2000 ) (see 
also their adaptively smoothed ROSAT image) who propose 
that rising radio bubbles (such as the X-ray holes) are re- 



Bruggcn & Kaiser 2001 



McNamara et al. 2001 



Quilis, Bower 



sponsible for the overall spiral structure. The outer surface 
brightness drop is traced by a rise in temperature (with in- 
creasing radius). This is reminiscent, although not a s large 



of the temperatur e drop at the cold front seen by Marke- 
vitch e al. (200C). If thermal conduction is an important 



effect in the intracluster medium, such fronts could reflect 
the magnetic field structure. 

A comparison of the metallicity distribution with the 
galaxy distribution from an optical image from the Digital 
Sky Survey does not show any obvious correlation between 
the galaxies and regions of higher metallicity. Projection ef- 
fects and the resolution of the Chandra metallicity map, 
however, do not allow to rule out any correlation yet. 

We find that the spectra extracted from the two central 
X-ray holes as well as the western X-ray hole are best-fit by 
hotter and more metal rich gas than their immediate sur- 
roundings. Using a spectral deprojection under the assump- 
tion of spherical symmetry in a sphere segment containing 
the northern inner X-ray hole, we have tested whether this 
could be due to a projection effect. Interestingly, we find that 
most of the X-ray emission in the hole can be explained by 
the projected emission of the shells further out. This directly 
addresses the issue of the X-ray gas content of the X-ray 
holes mentioned in the introduction; we find tight limits on 
the presence of an isothermal component in the X-ray hole, 
ruling out volume-filling X-ray gas with temperatures below 
11 keV at 3cr. 

The temperature distribution of the gas in the vicinity 
of the X-ray holes is of great importance for theories of ra- 
dio source heating to counter the short cooling time of the 
X-ray gas in this cluster (e.g., F00). Longer exposure times 
are needed to reveal if there is cool gas associated with the 



wake o f the bub ble, as, for example, in the model by Chura- 



zov et a l. (2000 ), who proposed that the radio bubbles would 
transport cooler gas from the cluster centre 'upwards'. The 
bubbles may also transport or drag metals upwards. The 
subject of radio bubbles and their interaction with the intr- 



F02). The Chandra 
results provide essential observational input to test these 
theories. 



ACKNOWLEDGMENTS 

We thank the anonymous referee for useful comments. ACF 
acknowledges support by the Royal Society. 



REFERENCES 

Allen S.W., Fabian A.C., Johnstone R.M., Arnaud K.A., 

Nulsen P.E.J., 2001, MNRAS, 322, 589 
Anders E., Grevesse N., 1989, Geochimica et Cosmochimica 

Acta, 53, 197 

Arnaud K.A., 1996, Astronomical Data Analysis Software 
and Systems V, eds. Jacoby G. and Barnes J., pl7, ASP 
Conf. Series volume 101 

Balucinska-Church M., McCammon D., 1992, ApJ, 400, 
699 



Blundell K., Kassim N, Perley R., 2000; preprint astro- 



ph/0004005 



Boehringer H., Voges W., Fabian A.C., Edge A.C., Neu- 
mann D.M., 1993, MNRAS, 264, L25 

Briggs S.A., Snijders M.A.J., Boksenberg A., 1982, Nature, 
300, 336 

Briiggen M., Kaiser C.R., 2001, MNRAS, 325, 676 

Churazov E., Forman W., Jones C., Bohringer H., 2000, 
A&A, 356, 788 

Churazov E., Briiggen M., Kaiser C.R., Bohringer H., For- 
man W., 2001, ApJ, 554, 261 

Conselice C.J, Gallagher J.S., Wyse, R.F.G., 2001, AJ, 
122, 2281 

de Young D.S, Roberts M.S., Saslaw W.C, 1973, ApJ, 185, 
809 

Dickey J.M., Lockman F.J., 1990, ARAA, 28, 215 
Fabian A.C., Hu E.M, Cowie L.L., Grindlay J, 1981, ApJ, 
248, 47 

Fabian AC, 1994, ARA&A, 32, 277 
Fabian AC, et al, 2000, MNRAS, 318, L65 (F00) 
Fabian AC, Celotti A., Blundell K.M., Kassim N.E, Per- 
ley R.A, 2002, MNRAS, 331, 369 (F02) 
Giovannini G., Cotton W.D., Feretti L., Lara L, Venturi 

T, 1998, ApJ, 493, 632 
Heinz S, Reynolds C.S., Begelman M.C, 1998, ApJ, 501, 
126 

Kaastra J.S., Mewe R, 1993, Legacy, 3, 16 

Liedahl D.A., Osterheld A.L., Goldstein W.H, 1995, ApJ, 

438, L115 
Lynds R., 1970, ApJL, 159, L151 
Markevitch M, et al, 2000, ApJ, 541, 542 
McNamara B.R, O'Connell R.W, Sarazin C.L, 1996, AJ, 

112, 91 

McNamara B.R, et al, 2001, ApJL, 562, L149 
Quilis V, Bower R.G, Balogh M.L, 2001, MNRAS, 328, 
1091 

Sanders J.S, Fabian AC, 2001, MNRAS, 325, 178 
Unger S.W, et al, 1990, MNRAS, 242, 33P 



10 Schmidt, Fabian and Sanders 

This paper has been typeset from a TffjX/ P^TfjX file prepared 
by the author. 



This figure "fig6_colour.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/0207290vl 



